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Abstract. - We study a qubit undergoing Landau-Zener transitions enabled by the coupling 
to a circuit-QED mode. Summing an infinite-order perturbation series, we determine the exact 
nonadiabatic transition probability for the qubit, being independent of the frequency of the 
QED mode. Possible applications are single-photon generation and the controllable creation of 
qubit-oscillator entanglement. 



Superconducting loops are promising candidates for solid state qubit implementations 
[1-4]. Since the direct observation of Rabi-oscillations in these systems [1], they form the 
basis of many experiments on coherent quantum dynamics. Particularly interesting are the 
experiments in circuit quantum electrodynamics (QED) [2-4], which is the solid-state analogue 
of a two-level atom in an optical cavity. Superconducting circuits possess the advantage that 
many of their parameters are tunable over a broad range. This can be exploited for controlling 
efficiently the qubits. 

One particular way of controlling a qubit works by switching the difference of its diabatic 
energies from a large negative to a large positive value, yielding an avoided crossing for the 
adiabatic energies. For sufficiently slow switching, the qubit will adiabatically follow its instan- 
taneous eigenstates. In the opposite limit of fast switching, however, the qubit will abandon 
the adiabatic eigenstate and undergo a so-called Landau-Zener (LZ) transition. The LZ tran- 
sition probability can be raised upon increasing the switching rate [5-7]. LZ transitions can 
be used to effectively control qubit gate operations [8,9] and to read out qubits [10]. Recently, 
LZ transitions have been observed in various experiments with superconducting qubits [11-14] 
and nanomagnets [15]. 

Landau-Zener transitions can also occur for a qubit that is coupled to a circuit oscillator. 
Then, the adiabatic following to the final ground state takes place even in the absence of a 
direct coupling between the two qubit levels, induced instead by the indirect coupling to the 
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oscillator. An oscillator in a highly excited coherent state can be described classically, yielding 
a non-monotonic behavior of the LZ transition probability versus the coupling strength [16]. 

Experimentally relevant is also the opposite limit in which the oscillator is initially in its 
ground state. Here, we shall investigate the Landau-Zener dynamics of such a qubit-oscillator 
setup and make specific suggestions for exploiting it to manipulate the quantum states of 
the qubit and the oscillator. We derive an exact analytical expression for the transition 
probabilities and we propose the creation of single circuit photons. 

In tunneling representation, a circuit QED setup is described by the Hamiltonian [2,3] 

H(t) = -ya, - + Mltfb + 7 (&t +b)[a x + l- 2N g ] (1) 

for the qubit, the circuit oscillator, and their mutual coupling. The cavity is modelled as a 
harmonic oscillator with the annihilation and creation operators b and b\ while o~ x and a z 
denote Pauli matrices. The electrostatic energy E c \ = 4E C [1 — 2N g ] is determined by the 
charging energy E c and the tunable gate charge N g . The tunable flux penetrating the 
superconducting loop will be used to drive the qubit. The flux controls the Josephson energy 
Ej(t) — £Jj >max cos[7r$(t)/$o] ; where <E>o is the flux quantum. The two-level approximation 
underlying the Hamiltonian JJJ is valid in the charge regime E c 3> Ej. In order to minimize 
decoherence, one typically operates the qubit at the optimal working point N g = h, so that 
E e i = [17]. Henceforth, we restrict ourselves to this optimal working point. The LZ dynamics 
can then be realised by switching the flux $(t) in such a way that Ej(t) = —vt, with v > 0. 
The duration of this linear sweep has to be long enough, so that transition probabilities 
have converged and the finite time interval can be extended to t = — oo . . . oo in calculations 
describing the dynamics. Since the energy splitting of the qubit can even vanish, it is essential 
to treat the qubit-oscillator coupling beyond rotating-wave approximation. The temperatures 
in circuit QED experiments [2] justify the assumption that both the qubit and the oscillator 
are initially in their ground states, i.e. |4 r (— oo)) = |t,0), where er 2 |T) = T}- 

Nonadiabatic transitions. - For E c \ = 0, the states ||) and ||) are eigenstates of the qubit 
Hamiltonian \vta z and any transition between these states can only result from the coupling 
to the oscillator. The central quantity of interest is the probability P^±(t) = 1 — Pf_>i-(i) 
that the qubit has flipped to the state \\). In the following, we derive an exact expression 
for Pf_>f(oo) = J2 n \(1i n \U(oo, — oo)||,0)| 2 where U(t,t') denotes the time-evolution operator 
and \n) an oscillator eigenstate. We start by a transformation to an interaction picture with 
respect to the uncoupled qubit and oscillator, Uo(t) = exp(— iClwbt) exp(— j^vt 2 a z ). This 
yields the interaction-picture Hamiltonian 

H(t) = 7 (6V m + b e' int ) exp ( - ^vt 2 a z ) a x . (2) 

A perturbation expansion of the probability amplitude A n = (f, n|£7(oo, — oo)|f, 0) results in 
the series A n = J2k=o(~' 1 /^) 2kank ^ with the 2/c-th order contribution a n k equal to 

/oo rt 2 h r t 2 p 2k , y . , 

dt 2k / dt 2k -i... dtiexp iV(OAA + — . 
J-oo J -co £—i lh ' 1 

~ (3) 
Since each H{t) flips the qubit exactly once, only even orders of 7 appear in A n . The coef- 
ficients = ±1, for I — 1, . . . , 2k, stem from the sign in the time-dependent phase of the 
creation and annihilation operators; C„fc({A}) = C n k(^2k, ■ ■ ■ , Ai) = (n\ ■ ■ ■ |0), where the dots 
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denote the combination of 2k creation and annihilation operators b and w that corresponds 
to the sequence A2/C, . . . , Ai. An important simplification of the A-summation results from the 
fact that C n k = whenever more annihilation than creation operators act on the oscillator 
ground state |0). Thus, we need to consider only those A-sequences that fulfill the relation 

1 

A* > 0, VI < 2k. (4) 

t'=i 

For the further evaluation, we substitute in Eq. the times tt by the time differences 
ti = t£+i — t[, for £ = 1, . . . , 2k — 1, where t = t<2k- Thus, we insert ti=t — Y^yZt T( -' > so that 
the integral in (0 becomes 

/qq 00 2/c 2 X • fc 2 X 

dt / drafe^.-.dnexp [iQ^A^i-^ T t') + ^f i J2{ 2Tn - 1 { t ~ T e')- T 2i-i} 



~°° J ° 1=1 t'=t " 1=1 l'=2l 

The ^-integration results in the delta function 

k 2k 



(5) 



27r<5 (iE^-i + fi E^)- ( 6 ) 

1=1 1=1 

From the inequality Q it follows that the second sum in the argument of the delta function is 
non- negative. Because the integration interval of all is [0 . . . 00), any non-zero contribution 
to the integral (0 comes from n = T3 = . . . = Tik-i = 0. Hence, the integral over the time 
differences T2, T4, . . . , T2k~2 must yield a distribution proportional to S(t{) 5(ts) ■ ■ ■ 5{t2U-i)- 
Evaluating the integrals over all T21 separately, one finds that such a distribution is obtained 
only if z2i'=i = for all £ = 1, . . . , k — 1. These fc — 1 relations together with the delta 
function © lead to the conditions X2i + Me-i = for all £ = 1, . . . , k. In combination 
with Eq. Q they imply that only those integrals with \g = — (— 1)^ are non- vanishing. In 
other words, we obtain the selection rule that to the final occupation probability only those 
processes contribute in which the oscillator jumps (repeatedly) from the state |0) to the state 
|1) and back. Hence the only relevant C n k reads C n k{— 1, 1, ■ ■ ■ , — 1, 1) = {n\(btf) k \0) — <5„,o- 
The remaining multiple integrations are performed as detailed in Refs. [18, 19] and result 
in a n k — 5 n fi(Trh/v) k /fc!, so that A n — S n .o exp(— ir^/hv). This implies that the oscillator 
returns to its ground state provided the qubit ends up in state If). Thereby, we arrive at a 
first central result, namely the exact transition probability 

P W (00) = 1 - P M (00) = 1 - e - 2 ^ 2 ' hv . (7) 

Most surprisingly, it does not involve the oscillator frequency f2. Because of this, Eq. (JJJ) 
allows a simple interpretation: if the Josephson energy is switched slowly (and the condition 
for that is Hv <C 7 2 ), then the qubit will follow the adiabatic ground state which at large 
times is the state ||). For large ftw/7 2 , the qubit will remain in state |f), corresponding to a 
nonadiabatic transition. Notice the important difference to the standard LZ problem that an 
intrinsic coupling between the qubit levels is replaced here by the qubit-oscillator coupling 7. 

Since only the states |f,0) and ||, 1) contribute to the perturbation series for Pj^f (00), it 
is tempting to conclude a posteriori that for the qubit-oscillator coupling, a rotating- wave ap- 
proximation (RWA) is justified. This would restrict the whole dynamics to the two mentioned 
states and indeed the same probability P|^|(oo) would be obtained. However, below we will 
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Fig. 1 - Landau-Zener dynamics for the coupling strength 7 = Q.&yfhv for various cavity frequencies 
fi. The dashed line marks the f2-independent, final probability J7J to which all curves converge. 



find that at intermediate times, states |T,2n) and |J,,2n+l) with n > can be considerably 
populated — the latter even at t — > 00. This demonstrates the surprising fact that the RWA 
yields the correct transition probability (JJJ even when the RWA is not justified. 

In order to gain information about the dynamics at intermediate times, we numerically 
integrated the Schrodinger equation for the Hamiltonian (|5J). The time evolution of the prob- 
ability that the qubit is in state j) is depicted in Fig.^ It demonstrates that at intermediate 
times, the dynamics depends strongly on the oscillator frequency f2, despite the fact that this 
is not the case for long times. For a large oscillator frequency, P\-,i (t) resembles the standard 
LZ transition with a time shift Ml/ v. 

If the system starts out not in the ground state but in a state IT,") with n > 0, some 
transition probabilities can still be obtained analytically. We find -Pj„_ »f( n + m ) = f° r ai l 
m > 0, being a special case of the selection rule derived in Refs. [20,21]. Moreover, P^ n ^ n = 
exp [— 27r(2n + l)j 2 /Hv] , indicating that |T, 0) is the most stable state. 

Manipulation of the oscillator state. - Having studied the qubit dynamics, we next focus 
on the modification of the oscillator state. Since we start out in the ground state |T,0) and 
the Hamiltonian (J2J correlates every creation or annihilation of a photon with a qubit flip, 
the resulting dynamics is restricted to the states |T,2ra) and |j,2n+l). Figure El reveals that 
the latter states survive for long times, while of the former states only |T,0) stays occupied, 
as it follows from the relation that A n oc 5 n> o, derived above. Thus, the final state exhibits a 
peculiar type of entanglement between the qubit and the oscillator, and can be written as 

|*(oo)> - a/i-P w M|T,0) + ^P T ^(oo) (cxU,1) + c 3 ||,3) + ...), (8) 

where |ci| 2 + |c3| 2 + ... = l. This implies that only odd- photon states are occupied if the qubit 
ends in |j), representing a highly nonclassical oscillator state. Qubit and oscillator end up 
fully entangled, in the sense that after tracing out the oscillator states, no coherence between 
the qubit states |T) and J.) is left. 

While P-|-_,^(oo) is determined by the ratio j 2 /hv, the coefficients C2 n +i depend also on 
the oscillator frequency. Experiments in circuit QED have been performed for 7 <C Kl in the 
coherent regime [2]. There, c\ ~ 1 to a very good approximation. Hence one can control via 
v the final state to be any superposition of |T,0) and ||,1). In particular, in the adiabatic 
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Fig. 2 - Population dynamics of individual qubit-oscillator states for a coupling strength 7 = 0.6V~hv 
and oscillator frequency £7 = 0.5^v/h. 

limit vft/7 2 <C 1, the final state becomes This has the important physical implication 

of the creation of exactly one photon in the cavity, triggered by a Landau-Zener transition. In 
an experiment, the photon will subsequently leak out of the cavity. By exploiting these two 
processes, we propose the following four-step LZ cycle for single-photon generation: The first 
step is single-photon generation in the cavity via the adiabatic LZ transition | T , 0) — > ||, 1), 
brought about by switching the Joscphson energy sufficiently slowly. Second, the photon is 
released from the cavity via the (controlled) cavity decay — > ||,0). In the third step, 
another individual photon is generated via the reverse LZ sweep |J.,0) — > |f, 1). Fourth and 
finally, a further photon decay completes the cycle. 

This scheme for repeated photon generation via Landau-Zener cycles makes use of two 
advantageous properties of circuit QED: first, the artificial atom is fixed at an antinode of the 
cavity, so that the atom-cavity coupling remains at a constant and high value. Second, qubits 
are highly tunable so that LZ sweeps can be made from minus to plus an "atomic" frequency, 
and back. The outlined scheme presents an alternative to the proposals for single-photon 
generation put forward in Refs. [22,23]. The main advantage of the present scheme is its 
robustness against parameter variations when operating in the regime 7 <C Ml. 

Experimental realisation. - In practice, the cavity frequency fl and the qubit-oscillator 
coupling 7 are determined by the design of the setup, while the Josephson energy can be 
switched at a controllable velocity v — ideally from Ej = — 00 to Ej = 00. In reality, 
however, Ej is bounded by Ej^ max which is determined by the critical current. The condition 
-Ej.max > HQ is required so that the qubit comes into resonance with the oscillator sometime 
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Fig. 3 - Probability of single-photon generation P\i,i) as a function of hv/j 2 , for LZ sweeps within 
the finite time interval [—T,T] with T > T m i n chosen such that vT = 3M/2. The initial state is 
|Tj0). Shown probabilities are averaged within the time interval "^Ml/v and '^hQ/v, whereby the 
small and fast oscillations that are typical for the tail of a LZ transition are averaged out. 



during the sweep. Moreover, inverting the flux through the superconducting loop requires a 
finite time 2T m j n , so that v cannot exceed i> max = £ , j,max/2T m ; n . In order to study under 
which conditions the finite initial and final times can be replaced by ±oo, we have numerically 
integrated the Schrodinger equation in a finite time interval [— T, T]. Results are presented 
in Fig. El The first three curves in Fig. |3 correspond to the experimentally relevant regime 
where 7 <C Ml. Here the essential dynamics is a LZ transition from \1, 0} to |, 1) around the 
time Ml/ v. Especially for the smallest coupling constants shown, the probability for single- 
photon generation practically equals P|^|(oo) given by Eq. (0. Thus we find that finite-time 
effects do not play a role as long as 7 <C Ml. Our predicted transition probabilities based 
on analytical results for infinite propagation time are therefore useful to describe the finite- 
time LZ sweeps. Figure also illustrates that the probability for single- photon production is 
highest in the adiabatic regime hv/-f 2 <C 1. Here the typical duration of a LZ transition is 
27/1; [16,24]. So in the regime of interest, the sought condition for a "practically infinite time 
interval" is vT — £j. max > + 27. For the unrealistically large qubit-oscillator coupling 
^/Ml = 0.5, reliable single-photon generation is less probable. This is so because (i) the LZ 
transition is incomplete within [— T, T]; (ii) more than two oscillator levels take part in the 
dynamics and more than one photon can be generated, as depicted in Fig. |2 and (iii) the 
approximation of the instantaneous ground state at t — —T by ||,0) is less accurate. 

For the setup of Refs. [2,3], a typical cavity frequency is fl = 2ir x 10 9 Hz. The sign of 
the initial Joscphson energy -Ej !m ax = ^ttH x 10 10 Hz can be inverted within T = 1/is so that 
« max = 2ttHx 10 16 s~ 2 . For the qubit-oscillator coupling strength we assume 7/27rfi = 3xl0 6 Hz. 
Notice that 7 <C hft, so that generation of more than one photon is negligible. Then, by 
choosing v — 0.05w max , one finds |ci| 2 > 1 — 10~ 7 and a superposition with equal weights 
for which Pj^| (00) = |ci| 2 P|^|(oo) = 0.5 is obtained. By choosing v = v max instead, a 
nonadiabatic transition occurs with the probability i-f^j(oo) ~ 0.97, while the remaining 
probability again corresponds to single-photon generation. As another extreme case, for the 
much slower sweep velocity v = 0.01u max , one is in the adiabatic limit in which a single photon 
is created with a probability of |ci| 2 P|^|(oo) = 0.97. Thus for the final state, one can obtain 
any desired superposition of the states ||, 0) and ]J,, 1} by choosing a proper value of v. 
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In summary, we have shown that the coupling of a qubit to a circuit cavity mode induces 
Landau-Zener transitions upon variation of the penetrating flux. For an oscillator initially in 
its ground state, we derived an exact closed expression for the transition probability which 
possesses an appealing form: it only depends on the interaction strength between the qubit 
and the oscillator and, in particular, it is independent of the actual oscillator frequency fi. 
Thus, the LZ dynamics can be manipulated via the velocity at which the qubit levels cross. For 
a circuit QED setup, we found that both nonadiabatic and adiabatic sweeps are feasible and 
can be exploited for quantum state preparation. Repeated adiabatically slow Landau-Zener 
sweeps allow the controlled and robust creation of single photons. 
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